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ABSTRACT 



Aims. Our analysis is aimed at characterizing the properties of the integrated spectrum of active galactic nuclei (AGNs) such as the 
ubiquity of the Fe Ka emission in AGNs and the dependence of the spectral parameters on the X-ray luminosity and redshift. 
Methods. We selected 2646 point sources from the 2XMM catalogue at high galactic latitude (\BII\ > 25 degrees) and with the sum of 
EPIC-PN and EPIC-MOS 0.2-12 keV counts greater than 1000. Redshifts were obtained for 916 sources from the Nasa's Extragalactic 
Database. We excluded sources classified as HII regions, groups/clusters, star-forming/starburst galaxies. The final sample consists of 
507 AGN. Individual source spectra have been summed in the observed frame to compute the integrated spectra in different redshift 
and luminosity bins over the range 0<z<5. Detailed analysis of these spectra has been performed using appropriately normalized 
background spectra and exposure time weighted response and ancillary files. 

Results. We find that the narrow Fe Ko- line at 6.4 keV is significantly detected up to z=l. The line equivalent width decreases with 
increasing X-ray luminosity in the 2-10 keV band ("Iwasawa-Taniguchi effect"). The anticorrelation is characterized by the relation 
log(EW Fe ) = (1.66 ± 0.09) + (-0.43 ± 0.07) log(L^_44), where EWf e is the rest frame equivalent width of the neutral iron Ko- line in 
eV and L x44 is the 2-10 keV X-ray luminosity in units of 10 44 erg s . The equivalent width is nearly independent of redshift up to 
z~0.8 with an average value of 101±40 (rms dispersion) eV in the luminosity range 43.5<logLje<44.5. Our analysis also confirmed 
the hardening of the spectral indices at low luminosities implying a dependence of obscuration on luminosity. 

Conclusions. We confirm that the neutral narrow Fe Ka line is an almost ubiquitous feature of AGNs. We find compelling evidence 
for the "Iwasawa-Taniguchi effect" over a redshift interval larger than probed in any previous study. We detect no evolution of the 
average rest frame equivalent width of the Fe Ka line with redshift. 

Key words, galaxies: active - X-rays: galaxies - quasars: emission lines 



1. Introduction 



have been significantly detected and characterized only in a few 



Active galactic nuclei (AGNs) are believed to be powered by ac- 
cretion ont o supermassive black holes (SMBH) (ISalpetedl 1964c 
iLvnden -Bell 1969; Rees 1984). Most of the accretion power 
is released in the innermost region arou nd the central SMBH 
(IShakura & Sunvaevlll973c iPringlel 1 1 98 lh . where the emerging 
spectrum is significantly affected by stro ng relativistic effects 
dPage & Thomdl974MCunninghamll975l) . The Fe Ka emission 
line is the most prominent spectral feature observed in the X-ray 
spectra of AGNs. The line profile serves as a powerful diagnostic 
tool to probe its origin and prov ides unique information on th e 
dynamics of the emitting region dFabian et al.ll 1 989c lLaorll 1 99 lh . 

The Fe Ko- emission was first identified as a common fea- 
ture in the X-r ay spectr a of Seyfert galaxies observed by Ginga 
(Pounds etal. 1990c iNandra & Pounds 1994). Early studies of 
type 1 AGNs with ASCA revealed resolved, broad and asymmet- 
ric Fe Ka lines in the majority of sources under study. The asym- 
metric li nes were described as being originated from the accre- 
tion disc (INandra etal .1119971) . Recent Chandra and XMM AGNs 
observations show that narrow Fe K emission is a ubiquitous 
feature in their X-ray s pectra (lYaqoob & Padmanabhanl 12004; 
iPounds & Reeves! l2002h . However, unambiguous broad lines 



Seyfe rt 1 galaxi es, examples being MCG-6-30- 15 (Fabian et al 
20021) . Mrk 205 dReeves et al.ll200ll) and Mrk 509 dPounds et al 
2001b . The number of spectral counts in AGNs is a major con 



straint in characterizing the line properties, since an accurate 
modeling of the underlying continuum is essential to measure 
the line parameters. 

Several studies have been conducted to charac terize the av- 
erage Fe Ko- emission in local and distant AGNs. INandra et al] 
(120071) performed a spectral analysis of a sample of 26 nearby 
Seyfert galaxies (z<0.05) observed by XMM-Newton. They 
found evidence for broad line emission in around 65% of the 
sample with typical equivalent wi dth (EW) being ~77±16 eV, 
when fitted with a broad Gaussian. iGuainazzi et al.l (l2006al) per- 
formed a similar analysis using a sample of 102 AGNs (includ- 
ing Seyfert galaxies, Radio Quiet Quasars) observed with XMM- 
Newton at z<0.5. They detected relativistic lines in about 25% 
of the sample objects (^50% when the "well exposed" spec- 
tra with > 10000 net counts in the 2-10 keV band are consid- 
ered). They also divided the sample in four equally populated 
luminosity classes and stacked residuals by modeling the con- 
tinuum in the 2.5-15 keV energy range with an absorbed power 
law and Compton reflection component. They found that the av- 
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erage EWs derived from the stacked spectral residuals analysis 
are <150 eV for all luminosity classes. They reported that the 
strongest relativistic profiles may be present in low-luminosity 
(L x <10 43 erg s" 1 ) AGNs. 

IStreblvanska et all (120051) derived an average rest frame 
spectrum of AGNs detected in a 770 ksec XMM-Wwtofi ob- 
servation of the Lockman Hole field. Using a sample of 104 
X-ray sources (53 type 1 and 41 type 2 AGNs and 10 galax- 
ies) covering 0<z<4.5 in redshift space, they reported a clear 
relativistic line in the average rest frame spectra with an EW 
of -560 e y and -460 e V for the type 1 and type 2 AGNs, re- 
spectively. ICorral et aD (120081) computed the mean Fe emission 
from a large sample of more than 600 type 1 AGNs spanning a 
redshift range up to -3.5. They detected significant unresolved 
Fe Ka emission line around 6.4 keV with an EW-90 eV and 
found no compelling evidence of any significant broad relativis- 
tic emission line in the average spectrum. 

iBrusa et alj d2005l) studied the average spectra of AGNs de- 
tected in the Chandra Deep Field North and South. They used 
a large sample of 352 spectroscopically identified AGNs and 
stacked X-ray counts in the observed frame in seven redshift 
bins over the range 0.5<z<4. The measured Fe Kg line EW s 
are consistent with those reported by IStreblvanska et al.l (120051) . 
However, a broad component was not required by the data. 
Based on these recent studies, the ubiquity of relativistically 
broadened Fe Ka emission in AGNs and redshift dependence 
of the line strength are still unclear. 

Studying the relation between emission line strength (char- 
acterized by the EW) and continuum luminosity is very impor- 
tant because it unveils the possible association b etween the emit- 
ting regions of these two spectral components. IB aid win! d 1977b 
first investigated the continuum luminosity dependence of the 
C IV ,11549 EW and discovered that the EW of the C IV line 
decreased with increasing UV luminosity, a phenomenon now 
known as the "Baldwin Effect". Since that time, confirming its 
existence and exploring its ori gin and evolution has become the 
subje c t of several studies (e.g . Dietrich et alj|2002t iGreen et all 



l200ll) . llwasawa & Taniguchil(ll993l) reported the anticorrelation 
between the Fe K emission line EW and X-ray luminosity, based 
on Ginga observations of AGNs. This relation was designated as 
"X-ray Baldwin effect" or "Iwasawa Taniguchi effect" (hereafter 
IT effect). This tre nd was confirmed by later ASCA observations 
of broad iron lines dNandra et al.1 19971) . In recent years, the exis- 
tence of the IT effect for the narrow Fe Ka emission lines in type 
1 AGNs obs erved by XMM-Newton and Chandra has been well 
established dPage et al.l200l IZhou & Wang||2005l: iBianchi et all 
l2007tlWu et al.l2009h . 

The most fundamental problem, however, is that the physi- 
cal causes for the "UV/optical Baldwin effect" and "IT effect" 
still remain unknown. Sever al models have been proposed to ac- 
count for these phenomena. Mushotzkv & Ferland ( 1984) sug- 
gested that the observed relation is due to a combination of 
decreasing ionization paramet er and covering fac tor with in- 
creasing continuum luminosity. iKorista et al.l (Il998l) pointed out 
that the Baldwin effect is driven by a softening of the ioniz- 
ing continuum and increasing metallicity towards higher lumi- 
nosities. Other studies propose that the underlying physical pa- 
rameters driving the Baldwin effect may be the Eddington ratio 



?: 



dBaskin & Laorll2004l: IWarner et al.ll2004tlDong et alll2009l) and 
black hole mass dXu et al.ll2008l) . lNandra et al.l (119971) attributed 



In the present study, we exploit the 2XMM catalogue to char- 
acterize the properties of the integrated spectrum to address the 
following key points : (i) ubiquity of the Fe Ka emission line 
in AGNs, (ii) origin of the IT effect and (iii) evolution of the Fe 
Ka equivalent width with redshift. The conventional approach 
to improve the counting statistics in AGNs involves unfolding 
the instrumental response by appropriately modeling the contin- 
uum and then shifting and adding the individual spectra, which 
in turn, maximizes the S/N ratio around the iron line energy. 
However, the main drawback is a model dependent parameter- 
ization of the underlying co ntinuum. We adopt an alternative ap- 
proach, first implemented bv lBrusa et alj d2005l) . in which source 
spectra are stacked in narrow redshift intervals. 

The paper is organized as follows. In the following Section, 
the selection criteria and properties of the sample are discussed. 
Section 3 outlines the spectral analysis. The following Sections 
then present and discuss the results. Thoughout this work a cos- 
mology with Q m =0.27, Q A =0.73 and H =70 km s" 1 Mpc" 1 is 
used. 



2. The Sample 

We use the sample derived from the 2XMM catalogue, the sec- 
ond comprehensive catalogue of serendipitous X-ray sources 
from the European Space Agency's (ESA) XMM-Newton ob- 
servatory. The catalogue contains 246897 X-ray source detec- 
tions which relate to 191870 unique X-ray sources dWatson et al.l 
120091) . We selected 2646 point source^] satisfying the follow- 
ing selection criteria : (i) high galactic latitude (\BII\ > 25 de- 
grees) and (ii) sum of the EPIC-PN and EPIC-MOS counts in 
0.2-12 keV greater than 1000. We cross-correlated the list with 
the Nasa's Extragalactic Database and obtained redshift infor- 
mation for 916 sources (-35% of the sample). 

Sources in the 2XMM catalogue are assigned 'automated' 
and 'manual' quality warning flags. These flags provide essential 
information about the possible problematic issues such as prox- 
imity to a bright source, a location within an extended source 
emission, insufficient detector coverage of the point spread func- 
tion of the detection etc. Additionally, the catalogue gives a sum- 
mary flag which combines the automated and manual quality 
warning flags. The summary flag, with values in the range 0-4, 
indicates the overall quality of each detection ordered by increas- 
ing severity. For a detailed description of possible values of the 
summary flag see lWatson et al.l d2009l) . 

The summary flags of the 916 sources were carefully 
checked. To obtain the cleanest possible sample, we retained 
the sources with '0' summary flags and removed the remain- 
ing sources from the sample. A source is assigned '0' summary 
flag when none of the automated or manual quality flags is set to 
True. We also used optical spectroscopic information for further 
screening of the sample and excluded all sources classified as 
HII regions, in groups/clusters and starburst/star forming galax- 
ies. The final sample consists of 507 point sources. 

We then grouped the sources according to their redshifts in 
13 redshift bins in the range 0<z<5. The choice of bin sizes is 
driven by a trade-off between the sum of net counts in each bin 
and to avoid too wide redshift bin preventing the detection or 
smearing the narrow line feature. The 0.2-12 keV flux (observed 
frame) given in the 2XMM catalogue was used to calculate the 



the IT effect to the presence of an ionized skin on the accretion 
disc, with the deg ree of ionization increasing with luminosity 
dNavakshinll2000allbT) . 



1 In the 2XMM catalogue point sources are characterized by the 
source extent parameter (EP_EXTENT=0). It is worthwhile to point out 
that the source extent below 6" is considered as a point source and the 
extent is re-set to zero. 
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Fig. 1. Luminosity-redshift distribution (top) and redshift distri- 
bution (bottom) of the sources in the final sample. 



Table 1. Statistics of the sources 



Redshift 


Total 


PN net counts 


Total 


bin 


sources 


(2-10 keV) 


Lbins 


0.0-0.2 


130 


78857 


6 


0.2-0.3 


42 


33896 


2 


0.3-0.4 


37 


15594 


1 


0.4-0.6 


47 


24349 


2 


0.6-0.8 


60 


35550 


2 


0.8-1.0 


33 


14502 


2 


1.0-1.3 


60 


30263 


2 


1.3-1.7 


52 


35735 


2 


1.7-2.1 


12 


5401 


1 


2.1-2.5 


15 


9352 


1 


2.5-3.0 


6 


5098 


1 


3.0-4.0 


4 


6717 


1 


4.0-5.0 


9 


12611 


1 



luminosity of the sources. In TableQ] we report for each redshift 
bin: the total number of sources, EPIC-PN counts in the 2-10 
keV rest frame and the total number of luminosity bins in which 
the sources were further splitted. The distribution in luminosity 
and redshift of the final sample is shown in Figure Q] 



3. Data Analysis 

For each object, source specific products such as the source and 
background spectra, images in the different energy bands and 
response files were retrieved from the XMM-Newton Science 
Archive (hereafter XSA). Keeping in mind the limitations gen- 
erally associated with the automatic detection process, several 
safety checks were carried out before using the archival prod- 
ucts. 



First, we examined the quality of the archival products by 
extracting source, background spectra, and response matrices 
from the observation data files (ODF) for a few sources using 
the XMM-Newton Science Analysis Software (SAS, ver. 7.1.0). 
The sources were arbitrarily selected with '0' quality summary 
flag and varying EPIC-PN counts spanning a range of ~300 to 
5000 in the 0.2-12 keV band. We cross-compared these ODF 
extracted spectra with the pipeline processed archival spectra re- 
trieved from the XSA. We found a very good agreement between 
these spectra for all sources under study. 

Second, we visually inspected the EPIC-PN images of the 
sources in the final sample to check if there are additional 
sources in the background annulus. In a few images, background 
annuli were found to have additional sources. We excluded all 
the archival products of these crowded field images from the 
stacking analysis. These exercises ensured us the goodness and 
reliable quality of archival products. Only EPIC-PN data has 
been used in our work. 

3.1. Spectral fitting 

Individual source spectra have been stacked in the observed 
frame in 13 redshift bins over the range 0<z<5. To avoid the 
contamination of the stacked spectrum by sources covering a too 
wide range of luminosities in a redshift bin, sources are further 
splitted into luminosity intervals. The choice of luminosity in- 
tervals associated to a redshift bin is driven by the number of 
sources in that redshift bin. For this reason, sources in all red- 
shift bins at z>1.7 are kept in only one luminosity interval. In 
all luminosity interval(s) of a redshift bin, the contribution of 
the single source counts to the total counts (0.2-12 keV) was 
checked before stacking spectra. The sources with more than 
40% contribution were removed from the luminosity interval. 
Appropriate background, response and ancillary files were cre- 
ated using the FTOOLS routines mathpha, addrmf and addarf, 
respectively. The response and ancillary files were weighted for 
the exposure time. We refer to these count-rate averaged X-ray 
spectra as the "real stacked spectra". 

The real stacked spectra were rebinned to have at least 20 
count s per bin and analyzed using XSPEC (ver. 12.4.0: lArnaudl 
119961) . In all redshift and their associated luminosity bins, the 
real stacked spectra were initially fitted with a single absorbed 
power law in the rest frame 2-10 keV band. In about half of 
the sample we found a significant excess above the power law 
continuum in the energy range from 6.4/( 1 +z mflA ) to 6.4/( 1 +z m ,„) 
keV, where z m ,„ and z max are the bin boundaries. To investigate 
the presence of a Fe emission line, we added a narrow redshifted 
Gaussian component (model wabs(pow+zgauss)) with the line 
width fixed at cr=0. 1 ke V and the redshift fixed at the median 
of the bin. We noticed that the fit improved in all luminosity 
bins, but with varying significance (see Table 0. An F-test has 
been performed to test the significance of the narrow redshifted 
Gaussian line component. Given the well known limitations of 
the F-tes t in determining the p recise significance of a spectral 
line (see iProtassov et al.ll2002l) . we ran simulations in order to 
quantify the Fe line significance in the real stacked spectra (see 
Subsection l3.3l l. 

In all redshift and luminosity bins at z<1.7, the peak energy 
of the Fe line is allowed to vary, whereas, at z>1.7, the centroid 
of the Fe line was fixed at 6.4 keV to constrain the fit. Whenever 
the addition of the line component was not significant (its sig- 
nificance <90% according to an F-test), an estimate of the upper 
limit to the rest frame EW of the neutral Fe emission line was 
inferred. Errors on the Fe line EW were calculated based on the 
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Table 2. Results of spectral fitting in different redshift bins. 



^median t±\r±(z) 

(keV) 
(1) (2) 



Binwidth Total Spectral 

(keV) sources counts 

(3) (4) (5) 



logLjf(2-10keV) T Eqw 

erg s _1 (eV) 

(6) (7) (8) 



X 2 /dof 

wabs(pow+zgauss) 

(9) 



;rMof 

wabs(pow) 
(10) 



F-prob 
(11) 



zbin: 0.0-0.2 



0.02 
0.03 
0.06 
0.09 
0.12 
0.16 



6.26 
6.22 
6.01 
5.89 

5.71 
5.51 



0.21 
0.28 
0.83 
0.81 
0.75 
0.67 



20(8) 
20(9) 
23(4) 
21(6) 
22(3) 
24(6) 



10440 

8421 

9433 

13986 

11932 

24645 



41.58!™, 
41.81!™* 

42.4333 

42.83!™' 
43.28!™' 

A.-X 04+0.03 



93+ 008 

"•"-0.07 

1 04 + ™ 6 

1-U ^-0.05 



1.06+ 



^+0.09 
-0.08 
i co+O.lO 
i - J -0.09 

1 20 + ™ 8 

1 - z ' u -0.08 
i Q7+0.05 
1 -"-0.04 



398 + |j 
455!« 

163 !s2 

105!^ 

75+ 3 ° 



491.9/435 
420.8/351 
373.4/385 
508.9/518 
474.1/472 
712.9/737 



642.5/437 


>99.999 


544.9/353 


>99.999 


399.3/387 


>99.999 


527.7/520 


99.991 


488.3/474 


99.906 


729.3/739 


99.976 













zbin: 


0.2-0.3 










0.25 
0.24 


5.12 
5.14 


0.32 
0.34 


21(2) 
21(2) 


6415 
27481 


43 14 + ™ 5 

^-'• 1 ^-0.31 

44 18 + ™ 3 


1 QQ+0.12 

1.88_ 007 

i 79+O.O3 
1 • ' -0.02 


201+no 

iul -109 
7Q+31 
' ^-30 


259.5/260 
728.8/708 


268.3/262 
746.2/710 


98.717 
99.976 












zbin: 


0.3-0.4 










0.34 


4.79 


0.34 


37(5) 


15594 


43 76 + ™ 3 


1 oo +0.05 
X - 8/ -0.03 


6 n+49 
u? -49 


439.6/477 


444.8/479 


94.026 












zbin: 


0.4-0.6 










0.47 
0.51 


4.37 
4.24 


0.54 
0.48 


23(2) 

24 


7099 
17250 


AT. tO+0.04 
"■•"-0.17 
44 4.1+0.04 


i 71+O.O8 
1 -' 1 -0.06 

i 74+0.04 

lr/ MJ.06 


144+95 
<98 


281.5/274 
475.3/485 


287.7/276 
479.4/487 


94.933 

87.255 












zbin: 


0.6-0.8 










0.68 
0.68 


3.80 
3.81 


0.40 
0.40 


30(4) 
30(3) 


9296 
26254 


43.88!™* 
44 67 +a02 


1 90+ 007 

i - 7u -0.07 

1 81+0-03 

00 -0.03 


163 +87 
qo+39 
yJ -19 


297.8/330 
602.3/561 


305.7/332 
616.9/563 


98.672 
99.881 












zbin: 


0.8-1.0 










0.89 
0.91 


3.39 
3.36 


0.33 
0.32 


17 

16 


7142 
7360 


44 n+o.03 

^■ l '-0.12 

44 69+™ 3 
tt.u? 010 


1 90 + ™ 9 

i - 7u -0.07 
1 89 +0-06 
1 ' 8 °-0.06 


<223 
<135 


254.1/262 
263.8/264 


258.4/264 
265.1/266 


88.928 
47.559 












zbin: 


1.0-1.3 










1.16 
1.18 


2.97 
2.94 


0.38 
0.36 


30 
30 


12833 
17430 


44 45 + - 02 

^■^ -0.08 

45 02 +a02 


2 06 +008 
z..ou 005 

1 8^+0-05 
1-OJ_0.04 


<183 
<51 


325.9/354 
423.0/417 


329.8/356 
423.0/418 


88.076 
5.595 












zbin: 


1.3-1.7 










1.44 
1.46 


2.63 
2.60 


0.38 
0.38 


27 
25 


12255 
23180 


445 .+0.04 
tt.Jl „ 7 

45 3i + ™ 2 


1 OQ+0.13 

1 89 +007 


<174 
<110 


358.9/352 
468.6/441 


361.5/354 
473.4/443 


71.753 
89.606 












zbin: 


1.7-2.1 










1.89 


2.22 


0.23 


12 


5401 


45 03+™ 4 


1 69 +017 

1 - u? -0.08 


<193 


199.1/200 


201.9/201 


90.608 












zbin: 


2.1-2.5 










2.22 


1.99 


0.18 


15 


9352 


45.38!™* 


1 Q9+0.14 
1 - 7Z '-0.12 


<48 


291.2/279 


291.2/280 


0.000 












zbin: 


2.5-3.0 










2.88 


1.65 


0.16 


6 


5098 


45 56 + ™ 4 


2 00 +a20 


<115 


177.5/190 


178.0/191 


50.559 












zbin: 


3.0-4.0 










3.31 


1.49 


0.25 


4 


6717 


46.18!™ 4 


1 66 +016 

1 - uu -0.14 


<106 


224.9/222 


225.6/223 


59.089 












zbin: 


4.0-5.0 










4.31 


1.21 


0.14 


9 


12611 


46 32+™ 3 


1 76 +013 


<83 


275.2/284 


276.5/285 


74.598 



Notes: Col. (1) : Median redshift of the bin. Col. (2) : Expected energy of the neutral Fe Ka line. Col. (3) : Observed frame bin width at the 6.4 
keV position, defined as AE= 6.4/1 +z min -6A/l+z max . Comparatively large bin widths in the redshift bin 0.0-0.2 are due to the variation in z,„,„ 
and z,„ ov introduced by splitting the sources in different luminosity bins. Col. (4) : Total number of stacked sources in different luminosity bins 
in each redshift bin (in parenthesis number of sources with significant Fe Ka line). Col. (5) : Total net counts in the rest frame 2-10 keV band. 
Col. (6) : Rest frame X-ray luminosity measured in the 2-10 keV band. Col. (7) : Power law photon index quoted over the 2-10 keV band; errors 
are reported at the 90% confidence level for one interesting parameter. Col. (8) : Rest frame equivalent width of the Fe Ka line. Col. (9) : x 1 and 
degrees of freedom using a single absorbed powerlaw and a redshifted narrow Gaussian line. Col. (10) : x 2 an d degrees of freedom using a single 
absorbed powerlaw. Col. (11) : The F-test probability of the line detection, for F-prob<90% the upper limits on the EW are quoted at the 90% 
confidence level. 
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error in the line normalization considering the 90% confidence 
interval. The fitting results are summarized in Table|2] 

The first part of Table [2] lists for each redshift bin and its 
associated luminosity bins, the expected positions of the Fe Ka 
line, the bin width, the source statistics (number of sources and 
number of counts) and the average 2-10 keV luminosity. In the 
second part of Table [2] results of the spectral fitting in each red- 
shift bin (splitted into different luminosity bins) are presented, 
namely the photon index of the underlying power-law contin- 
uum (T), the EW of the Gaussian line component,^ 2 and degrees 
of freedom associated with an absorbed powerlaw and with an 
additional redshifted narrow line component. The F-test proba- 
bility of the line detection is also quoted. 

3.2. Contribution from single sources 

In redshift and luminosity bins where a significant Fe Ka line 
is detected in the real stacked spectrum, we performed spectral 
analysis of each source spectrum to check for the presence of a 
significant Fe Ka line. We first grouped each source spectrum 
to have minimum 20 counts per bin. As our main interest is to 
identify spectra with significant Fe Ka line, we restricted our 
spectral analysis in the 2-10 keV rest frame energy range. For 
simplicity we adopt the same continuum model for all spectra as 
for the real stacked spectrum. We then added a narrow (cr=0.1 
keV) redshifted Gaussian line with the redshift fixed for each 
source. We performed an F-test to quantify the significance of 
the Fe Ka line in each spectrum. We considered significant Fe 
Ka detection in the spectrum if the F-prob of the line detection 
was >90%. The number of sources found to have significant Fe 
Ka line are quoted in parenthesis in Table [2] 

We then investigated the Fe Ka contribution from spectra 
with and without significant Fe Ka lines by creating stacked 
spectra of both subsamples. These stacked spectra were analyzed 
with the same models as that used for the real stacked spectra 
and fit parameters were derived. The stacked spectra of sources 
without significant Fe Ka line detection show the presence of 
Fe Ka line with varying significance depending on the redshift 
and luminosity bin. Typical Fe Ka EWs of these spectra span 
a range from ~200 eV in the lowest redshift bin to ~100 eV in 
the redshift bin 0.6-0.8. In the top panels of Figures [2] and [3] we 
show spectral fits with respect to a single absorbed power law to 
three stacked spectra, from top to bottom, sources with signifi- 
cant Fe Ka line detection, real stacked spectra and sources with 
undetected Fe Ka line. In the bottom panels Data/Model ratios 
are shown for the real stacked spectra. 
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Fig. 4. Rest frame equivalent width as a function of redshift. 

Each simulated stacked spectrum was subsequently fitted 
with a single absorbed powerlaw and a redshifted Gaussian line 
with the line width and redshift fixed at 0. 1 ke V and Zmediam re- 
spectively. Then, histograms of all the spectral fit parameters 
were created and fitted with a Gaussian function to derive the 
mean and standard deviations of their distributions. Comparison 
of the spectral parameters of the real (see Tabled and simulated 
stacked spectra yielded a good agreement between them, in all 
cases within the statistical errors. Thus, we quantified that the 
EW of the Fe Ka derived from the spectral analysis of the real 
stacked spectrum represents the average of the EWs of the Fe 
Ka lines of the individual spectra comprising the real stacked 
spectrum. 

4. Results and Discussion 

Here we discuss the results of our work in different redshift and 
luminosity bins. At low redshifts (z<0.4) F-test probability of the 
Fe Ka line detection is greater than 90% in all luminosity bins. 
In the redshift range 0.4<z<0.8 the average F-test probability 
of the Fe Ka line detection is ~95% in the lowest luminosity 
bin. An upper limit on the Fe Ka EW was derived for the non- 
significant line detections in all redshift and luminosity bins at 
z>0.8. 

One key advantage of our sample with respect to previous 
published works is that it includes higher statistics covering a 
wide range of luminosities at various redshifts. We can therefore 
address the separate redshift and luminosity dependences of the 
Fe line strength. 



3.3. Simulations 

We performed extensive simulations to assess the accuracy of 
the stacking procedure. We created "simulated stacked spectra" 
in redshift and luminosity bins where a significant Fe Ka line 
is detected in the real stacked spectrum. To create the simu- 
lated stacked spectrum we considered separately the contribu- 
tion from sources with and without significant Fe Ka line. For 
the sources with significant Fe Ka line, the best fit parameters 
retrieved from the individual spectral analysis were given as in- 
put parameters. For the sources without significant Fe Ka line, 
the average values measured from their stacked spectrum (blue 
spectrum in Figures[2]and[3]l were used as input parameters. We 
produced 1000 simulated stacked spectra in each of the redshift 
and luminosity bins. Identical response and effective area files 
were used for the real and simulated stacked spectra. 



4. 1 . Evolution of the Fe line equivalent width with redshift 

First, we searched for a dependence of the Fe Ka line intensity 
on cosmic time (redshift). We selected the narrow luminosity 
range of 43.5<logLx<44.5 to cover a redshift range as large as 
possible for an almost constant luminosity. This sample contains 
206 sources grouped in 7 redshift bins in the range 0<z<1.3 and 
an average luminosity of logLx(2-10 keV)=44.04±0.30 (rms 
dispersion). We find that the intensity of the Fe Ka line does 
not change significantly up to z~0.8 as displayed in Figure [4] 
where the measured EWs with their associated errors and upper 
limits are plotted. The upper limits at z>0.8 are also consistent 
with a non evolution. The non evolution of the Fe Ka line in- 
tensity can be inte rpreted as constant Fe abundance with redshift 
dBrusa et al.ll2005l) . Studies of quasar elemental abundances (de- 
rived from broad optical emission lines) also reported a similar 
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i non evolutionary) tre nd of metalli city with redsh i ft up to z~4.5 
HamannetaT1l2007h . In contrast, iBalestra et all (120071) found 
significant evidence of a decrease in the average Fe abundance 
of the intra-cluster medium as a function of redshift, at least up 
to z^0.5. The comparison of these two results implies a different 
evolution of the Fe abundance with redshift in these two very 
different systems (AGNs and ICM) and/or very different scales 
(pc versus Mpc). 

4.2. The IT effect 

We then investigated for a dependence of the Fe Ka line strength 
on the X-ray luminosity. Figure[5]shows the relation between the 
average rest frame EW and 2-10 keV luminosity in all redshift 
and their associated luminosity bins. The weakening of the line 
strength with increasing luminosity is clearly present. To derive 
the relation between the tw o parameters we followed the method 
of dGuainazzi et al.ll2006bl) . Their method can be summarized as 
follows. First, a set of Monte-Carlo simulated data is derived 
from the EWs (both upper limits and measurements with their 
associated symmetric errors) according to the following rules: 
(i) each measured EW was substituted by a random Gaussian 
distribution, whose mean is the best fit value and whose stan- 
dard deviation is its statistical uncertainty; (ii) each upper limit 
(U) was substituted by a random uniform distribution in the in- 
terval [0,U]. Then, ordinary least squares fits were performed on 
each Monte-Carlo simulated data set. The mean of the slopes 
and the intercepts derived from the fits of each data set are taken 
as the "best fit parameters of the linear relation". We obtained 
the following best fit relation: 



log(EW Fe )=(1.85±0.1 1)+(-0.32±0.07) log(L x ,44) 



log(EW Fe ) = (1.66 ± 0.09) + (-0.43 ± 0.07) log(L x ,44) 



(1) 



where EWf e is the equivalent width of the narrow Fe Ka line 
in eV and Lx,44 is the 2-10 keV X-ray luminosity in units of 
10 44 erg s -1 . To test the significance of this anti-correlation, we 
calculated the Spearman's rank coefficient (Rs) for each Monte 
Carlo simulated data set. The mean value is taken as the R$ of 
Equation Q] The R5 for this correlation is -0.74, corresponding 
to a Null Hypothesis Probability (NHP) of 2.95xl0~ 4 . 

Our data therefore strongly imply that the Fe Ka line 
strength in the selected samples decreases as a function of in- 
creasing luminosity in all redshift bins confirming the IT effect. 
In the previous Subsection 14. 11 we explored the redshift depen- 
dence of the Fe Ka line intensity. The observed non evolutionary 
behavior assured us that no additional redshift effects are intro- 
duced in the IT relation given by EquationQ] which includes dif- 
ferent luminosities and therefore on average different redshifts. 
This indicates that the driver of the observed IT correlation is 
luminosity and not redshift. 

However, the slope of the EW Fe versus Lx relation we ob- 
tained is st eeper than the value reporte d in previou s studies (- 
0.20+ .0.03 (llwasawa & Taniguchill993h : -0.17+0.08 dPage et all 
12004 : -0. 17±0.03 dBianchi et al.ll2007h ). The normalization of 
our EWf,, versus Lx relation is statistically consistent with that 
reported bv lBianchi et al.l (l2007h . To check whether this discrep- 
ancy in the slope is due to the 2 extreme points at Lx<10 42 erg 
s -1 , we removed these points and subsequently performed the 
linear fits. We still obtained a slope of -0.43+0.09 with Rs=- 
0.67 and NHP=3.24xlO~ 3 . 

This discrepancy in the slope can possibly be attributed to 
th e difference in composition of the samples. Indeed, the sample 
in iBianchi et al.l (120071) consists of 157 local unobscured radio 
quiet AGNs (almost 90% within z<l), whereas our sample com- 
prises a mixture of 507 obscured and unobscured sources over a 
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Fig. 6. Rest frame equivalent width of the narrow Fe Ka line in 
the redshift interval 0<z<0.8 plotted against the X-ray luminos- 
ity in the 2-10 keV band. (Symbols, line styles and colors same 
as in Figure |5]). 



broader z range (69% at z<l and 31%atz>l). The difference in 
the slope we obtain for the whole sample can also be driven by 
the upp er limits in th e hig hest redshift bins which are not sam- 
pled bv lBianchi et alJ (120071) . When all bins in the redshift range 
0<z<0.8 are considered (see Figure |6), the inverse correlation 
between the EW and X-ray luminosity becomes 



log(EW Fe ) = (1. 85 + 0.11) + (-0.32 + 0.07) log(L x , 44 ) 



(2) 



with R s =-0.64 and NHP=3.63xl0~ 2 . The slope of the anticor- 
relation is statistically consistent in Equation Q] and [2] However, 
the best fit value is flatter in Equation [2] This in turn supports 
the hypothesis that upper limits are an important factor in the 
calculation of the fit relation. 

The physical origin of the IT effect is related to where the 
line is originated. Several studies have confirmed that the narrow 
neutral Fe Ka lin e is a ubiquitous component in AGNs spectra 
dPage et al.l l2004t lYaqoob & Padmanabhan] l2004t ICorral et al.l 
120081) . The measured widths vary substantially from source 
to s ource and cover a wide range from 1000 to 15,000 km 
s" 1 (lYaqoob & Padmanabhan! 12004 . thus being consistent ei- 
ther with an origin from the broad line region (BLR) or the 
torus (obscuring circumnuclear matter) envisa ged in orientation- 
depen dent unification sche me for AGNs (lUrrv & Padovanj 
1 19951) . The original IT effect dlwasawa & Taniguchill 19931) was 
attributed to a decrease in the covering factor of the BLR with 
increasing X-ray luminosity. However, recent studies report that 
no strong correlation exists between t he Fe Ka cor e width and 
prominent BLR H/3 ,14861 line width dNandrall2006l) . There is a 
similar lack of co rrelation betwee n the EWs of the Fe Ka and 
CIV ,11549 lines dWu et al.ll2009l This implies that the Fe Ka 
line is unlikely to be produced in the BLR. 

Based on previous studies, the abundance effect and black 
hole mass can be ruled out as the primary physical para meters 
driving the correlation (Page et al. 2004; Zhou & Wang 2005; 
lNandrall2006l: IBianchi et alJl2007l) . The absence of evolution of 
the Fe Ka EW with redshift (see Subsection 14. U further con- 
firms this hypothesis. Assuming that the line is originated in the 
torus, the luminosity dependent covering factor of the torus turns 
out to be the most likely e xplanation for the IT effect (Page et al. 
12004 IBianchi et alj2007l) . A lower covering factor signifies that 
the torus obscures a smaller solid angle around the nucleus, and 
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Fig. 7. The variation of the powerlaw photon index as a function 
of the rest frame 2-10 keV luminosity is illustrated. Symbols 
reflect different redshift bins as in Figure [5] 



vious studies, likely due to the difference in composition of 
the samples. The exact physical cause of the IT effect is un- 
known, one promising explanation is a decline in the cover- 
ing factor of the putative molecular torus as the luminosity 
increases. 
• We find a hardening of the spectral indices at low luminosi- 
ties in the redshift bin 0.0-0.2. This can be linked to the de- 
pendence of obscuration with luminosity. 

Large area X-ray surveys from sensitive next generation 
satellites (e.g., eROSITA and WFXT), when coupled with deep 
multiwavelength photometry (e.g., PanSTARRS and LSST) and 
spectroscopic campaigns (SDSS, LAMOST, and Euclid) will al- 
low to study the dependence of the EW strength with redshift and 
luminosities on samples of several order of magnitudes larger 
number statistics providing a unique sampling of the L-z plane. 



hence less Fe K-edge photons are captured and generate less Fe 
Ka photons, leading to a reduction in the EW. In this way, an 
increase in luminosity and a subsequent decrease in the covering 
factor could explain the observed anticorrelation of the EW and 
X-ray luminosity. 

An immediate implication of the luminosity dependent cov- 
ering factor of the torus is the variation in the underlying con- 
tinuum with luminosity as high covering factor at low luminosi- 
ties would cause spectral hardening by suppressing the soft part 
of the spectrum. To confirm this argument, we investigated the 
dependence of the power law photon index (characterizing the 
underlying continuum) on the X-ray luminosity in different red- 
shift bins. In order to avoid the complexity of the AGNs spectra 
at energies below 2 keV, we performed spectral analysis in the 
2-10 keV rest frame band. Figure [7] illustrates the variation of 
the power law photon index (F) with the 2-10 keV luminosity in 
different redshift intervals. We found that the T in our stacked 
spectra sample is hard at low X-ray luminosities (Lx<10 43 erg 
s ) sampled mostly in the redshift bin 0.0-0.2 (see also Figure 
2). A hard T is indeed expected in the case of increasing absorp- 
tion at these luminosities due to a high covering factor. The trend 
also reflects the decrease of the fraction of the absorbed sources 
with the X-ray luminos ity revealed from extensive works in hard 
X-ray surveys studies dUeda et al.ll2.Q03b lLa Franca et al.l l2005: 
lHasingerll2008l) . 

5. Conclusions 

We have compiled a large sample of 507 sources selected from 
the 2XMM catalogue covering the redshift range 0<z<5 and 
nearly 6 orders of magnitude in the X-ray luminosity. To im- 
prove the signal-to-noise ratio and investigate the average X-ray 
spectral properties of the sources, we computed the integrated 
spectra in narrow intervals of redshift and luminosity. We per- 
formed simulations to assess the accuracy of the stacking proce- 
dure. The main results can be summarized as follows. 

• We confirm that the narrow Fe Ka line is an almost ubiqui- 
tous feature in the X-ray spectra of AGNs. 

• In the integrated spectra spanning redshift range 0<z<1.3 at 
nearly constant luminosity we detect no strong trend in the 
Fe Ka line EW with redshift. 

• We find a significant anticorrelation between the Fe Ka EW 
and 2-10 keV rest frame luminosity across a wide redshift 
range not probed in prior studies. The slope of the anticorre- 
lation we obtained is steeper than the value reported in pre- 
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Fig. 2. Top panels in each figure show spectral fits (mo wabs(pow)) to three stacked spectra in the quoted redshift bin, from top 
to bottom, sources with significant Fe Ka line detection, real stacked spectra and sources with undetected Fe Ka line. For visual 
clarity, the stacked spectra of sources with and without significant Fe Ka line detection have been rescaled by the factor listed in 
parenthesis. In the bottom panels Data/Model ratios are shown for the real stacked spectra. The vertical line is drawn at the expected 
position for the redshifted 6.4 keV Fe ka line, while the shaded region encompasses the bin width reported in Table 2. The average 
rest frame luminosities (2-10 keV) obtained for the best fit model (a single absorbed power law and a narrow redshifted Gaussian 
line) are also quoted for the real stacked spectra. 
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Fig. 3. Representation of the data in this figure is the same as in Figure [2] 
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Fig. 5. Rest frame equivalent width of the narrow Fe Ka line (<x=0. 1 keV) in different redshift intervals plotted against the X-ray 
luminosity in the 2-10 keV band. The solid black line indicates the best fit to all the detections including u pper limits. The fit relation 
is reported on the top. The solid green line represents the best fit relation reported bv lBianchi et alj d2007l) . The combined error range 
on the slope and normalization of the best fit is also shown by dotted lines in the respective colors. 



